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Quantum size effects in armchair graphene nano-ribbons (AGNR) with hydrogen termination are 
investigated via density functional theory (DFT) in Kohn-Sham formulation. "Selection rules" will 
be formulated, that allow to extract (approximately) the electronic structure of the AGNR bands 
starting from the four graphene dispersion sheets. In analogy with the case of carbon nanotubes, a 
threefold periodicity of the excitation gap with the ribbon width (N, number of carbon atoms per 
carbon slice) is predicted that is confirmed by ab initio results. While traditionally such a periodicity 
would be observed in electronic response experiments, the DFT analysis presented here shows that 
it can also be seen in the ribbon geometry: the length of a ribbon with L slices approaches the 
limiting value for a very large width 1 <C N (keeping the aspect ratio small N <^ L) with 1/N- 
oscillations that display the electronic selection rules. The oscillation amplitude is so strong, that 
the asymptotic behavior is non-monotonous, i.e., wider ribbons exhibit a stronger elongation than 
more narrow ones. 

PACS numbers: 73.21. Hb,73.22.Pr,61.48.De 

Keywords: Graphene, Nano-ribbon, Quantum Size Effects, Edge Effects 



The nonstandard electronic properties 0, Q along with 
improved fabrication techniques have moved graphene 
and its allotropes into the focus of frontier research in 
recent times. [3| The presence of an edge in graphene 
nano-ribbons (GNR) takes a significant influence on this 
electronic structure. 0-0 In contrast to pure graphene, 
GNRs with a proper edge exhibit a finite bandgap poten- 
tially useful for device applications. Therefore, graphene 
nano structures with well defined orientation and edges 
have become a research field of their own. 7 19] 

As is well known, zigzag-edged ribbons have two flat 
metallic bands near the Fermi energy possibly lead- 
ing to magnetism. By contrast armchair ribbons 
(AGNR) are semiconducting with width-dependent band 
gap s and without affinity to magnetic instabilities. [il. [23l— 
|31| In this letter we further investigate the electronic and 
atomistic structure of mono-hydrogenated AGNRs and 
show that interesting quantum effects arise, nevertheless. 
We formulate selection rules for the transverse momenta 
of the ribbon. They identify those lines in the (extended) 
Brioullin-zone of graphene that resemble the electronic 
structure of AGNRs of a given width N, the number of 
carbon atoms in transverse direction, see. Fig. [T] In this 
way we infer that a reasonable first approximation for all 
energy bands of an AGNR is encoded in a single selection 
rule. Similar to the case of a carbon-nanotube the selec- 
tion rule predicts a three-fold periodicity of the bandgap 
in N. Because the same selection rule can be applied to 
all bands, one might suspect this periodicity to appear 
also in the atomic structure of the ribbon. Indeed, this is 
what has been observed previousl y in the edge stress and 
energy for nonpassivated ribbons. |32| Our detailed DFT- 
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FIG. 1. (Color online) Left: Atomistic structure of armchair 
graphene nanoribbons (AGNR). They come in two species 
with odd and even number of carbon atoms N in transverse 
direction. Right: The Brioulline zone of graphene appropri- 
ately parameterized for locating energy bands of the AGNR. 



anatomy of AGNRs reveals however that the three-fold 
periodicity also appears in the atomic structure of hy- 
drogen terminated AGNR: the longitudinal deformation 
of the unit cell of an AGNR with respect to the bulk 
graphene value is described in leading harmonic approx- 
imation by a term ~ cos(2irN/3)/(N — 1). 

Our results imply that quantum size effects can be 
studied experimentally in AGNR by atomic structure 
determination, namely by comparing the length of AG- 
NRs (having the the same number of carbon slices, L) of 
neighboring width N, N + 1, . . .. This is in marked con- 
trast to traditional approaches in meso- and nanoscopic 
structures that investigate quantum size effects near the 
band-gap by directly probing the electronic excitation 
spectrum, e.g. in transport measurements. 
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FIG. 2. (Color online) Bandgap A(TV) of hydrogen terminated 
AGNRs as a function of 1/N. 



Method and model. All DFT calculations pre- 
sented in this letter are performed using a plane-wave 
basis set and the projector-augmented wave (PAW) 



method. [331. |34( as implemented in the Vienna Ab ini- 
tio Simulation Package (VASP).[35T437l| Based on a com- 
parison of the graphene lattice-parameter obtained using 
several local, semi-local and meta-semi-local functionals 
(details are available in the supplemental material), the 
exchange-correlation functional of Perdew and Wang was 
chosen for the calculations. (3^ For all geometries, the 
atomic structure is fully optimized using a conjugate gra- 
dient algorithm. A convergency criterium of 1 meV/A is 
used for the forces on the atoms. [39J 

Electronic structure of H terminated AGNR. In 
Fig. [5] the gap for electronic excitations of an AGNR is 
plotted over the inverse system width. The data splits 
into three sets exhibiting an oscillatory behavior. [3(J For 
TV+1 divisible by three the smallest gaps A(TV) are seen 
and interpolate smoothly with 1/N into the bulk limit 
A=0. For N divisible by three the gaps are larger, also 
interpolating smoothly into A=0. Largest gaps are en- 
countered in the remaining case, N— 1 divisible by three. 

This overall phenomenology is reminiscent of the situ- 
ation with carbon nanotubes (CNT)[4tJ and we will in- 
vestigate it now closely. Zig-zag CNTs are similar to 
AGNR's with the H-termination replaced by periodic 
boundary conditions. The tube's electronic structure is 
understood by imposing selection rules for allowed trans- 
verse wavenumbcrs on graphene's band-structure in re- 
ciprocal k-space. These selection rules reflect the cylin- 
drical geometry. Three classes of armchair CNTs are thus 
obtained. (40l| In a very similar way, also the electronic 
structure of AGNR can be understood as a cut through 
graphene's dispersion sheets along lines of selected trans- 
verse momenta k y . 

To find the appropriate selection rule for the AGNR's 
we map the (valence) bands e(k x ) of the N = 5 to 9 rib- 
bons onto the four (valence) dispersion sheets, e(k x ,k y ), 
of pure graphene. The sheets are plotted as color maps in 




FIG. 3. (Color online) Energy bands for the N = 6 AGNR 
(long dashed green lines) plotted into the graphene dispersion 
sheets in the \t x , k H plane. The dashed white lines indicate 
estimates based on the selection rule, see Eq. JTJ. 



Fig. [3J For each sheet, k y values in graphene's Brioullin 
zone were selected such that the obtained set of €k y (k x ) 
functions most closely resembles the band-structure of 
the AGNR. (The comparison of the functions for fitted 
k y to the actual ribbons bands is shown for all four sheets 
of the N = 5 to TV = 9 AGNRs in Figs. 1 to 10 of the 
supplementary material.) The fitting process produces 
TV bands in sheets 1, 2 and 4. Due to the H-atoms at 
the edges the third sheet accommodates two more bands, 
N + 2. 

We infer from Fig. [3] that in sheets 1, 2 and 4 the 
AGNR-bands are approximately equidistant with Ak y sa 
fcmax/ (TV+1) which motivates the following selection rule: 
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where i runs from 1 to TV for the sheets 1 to 3 and from 
to TV + 1 for sheet 4 to also include the two £f-bands. 

The AGNR-bands as predicted by the rules, Eq. (P) 
are also indicated in Fig. [3j the root mean square devi- 
ation per sheet of the fitted values to the selection rule 
predictions are detailed in Fig. [4] We find a good match 
for the uppermost sheet 4, especially in the vicinity of 
the K-point. The overall discrepancy between the pre- 
dicted and the true transverse momenta of the ribbon 
bands never exceeds 15% of the interband spacing, ex- 
cept for the smallest ribbon widths. An exception to this 
rule is presented by the third sheet. Here, the presence of 
the 7i-atoms interferes and the selection rule gives only 
semi-quantitative information. 

In analogy to the case of carbon nanotubes, also the 
proposed rules ((T|) allow for a qualitative understanding 
of the electronic structure. Indeed, for the AGNR with 
mod (TV + 1, 3) = the sheets are divided into a multi- 
ple of three equally wide sections. Hence, the k y = 1/3 
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FIG. 4. (Color online) Root mean square deviation per sheet 
of the fitted k y values to the values predicted by the selection 
rule. The deviations in sheets 1, 2, and 4 remain smaller than 
15% of the line spacings independent of the ribbon width 
N. The bigger deviations in the fourth sheet result from the 
presence of the hydrogen edge. 
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FIG. 5. (Color online) Evolution of the edge C-C bonds plot- 
ted against 1/N. For the location of the bonds, al, a2, bl, 
and b2, see Fig. [T] The growth of the sample width is in 
negative y-direction, so a2 and b2 oscillate, while al and bl 
do not. The lower trace corresponds to the outer bonds, the 
upper trace to the (more) inner bonds. 



line goes though the K-point (fourth sheet, Fig. [3]) pre- 
dicting a ribbon with an anomalously small bandgap.[3fjj| 
After applying the selection rule to the other ribbons as 
well, we recover the values for the band gaps obtained by 
earlier tight binging calculations, with nearest neighbor 
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42j In addition the selection also works 
for all the other sheets. The particularities of these sheets 
is not only encoded in the K-point behavior but also in 
the position of Van-Hove singularities. Since the posi- 
tioning of the electronic levels with respect to these fea- 
tures also oscillates with period three, one expects that 
the oscillation behavior of the electronic structure also 
carries over to the atomic geometry. 

Atomic Structure of H Terminated AGNR Wc 
now investigate how the electronic structure translates 
into the atomistics of the AGNR, which we have also de- 
termined within our ab initio approach. For the cell ge- 
ometry we adopt the nomenclature as depicted in Fig. [1] 
The extension of the unit cell in the x, y-dircctions is in- 
dicated by d and w, respectively. For ease of comparison 
to bulk values we introduce the following dimcnsionlcss 
quantities: Ad = —A 1, Aw = , , , . w — jryr — 1, and 

H V3a ' (JV-l)(a /2) ' 

AA = (Ad + l)(Aw + 1) — 1, with ao the lattice parame- 
ter of bulk graphene calculated using the same functional 
and accuracy. All three become zero in the bulk limit, 
N -> oo. 

In the mono-hydrogen termination AGNR (H-AGNR) 
each carbon edge atom binds two other carbon atoms and 
one hydrogen atom. Considering the atomic structure of 
H-AGNR's this leads to a two-fold periodicity associated 
with TV even and N odd, see Fig. [1] AGNRs of con- 
secutive widths, N = 2n and iV = 2n+l, have different 
structural patterns; 2n— AGNRs exhibit re— 1 hexagons 
in all lateral sections while (2re+l)— AGNRs have a se- 
ries of alternating n and re— 1 hexagons. The effect of 



this pattern can be visualized, by plotting the edge C- 
C bond- lengths, for the fully relaxed structures over the 
reciprocal of the ribbon width 1/N, Fig. [5j The way the 
edge exerts a pressure on the ribbon is, however nontriv- 
ial. Although the overall effect is elongating, we observe 
that the outer C-C bonds al are always shorter than the 
inner ones bl (Fig. [l}. Even more, the sum of al and 
a2 is shorter than two times the optimal bulk C-C. The 
difference survives in the large N limit. The main cause 
of the elongation is hence not a change in bond lengths 
but the opening of the angles a and /3. 

The effect is also displayed in Fig. [6] which now empha- 
sizes healing of the geometry of the unit cell and allows 
for a simplified quantitative analysis: 

Ad = K d (N)/(N-l). 

The pre-factor is found to match the empirical form, 

Kd (N) ~ 4 0) (l + 4 15 cos (27TJV/3)) , (2) 

with fitting parameters that can be read of the inset of 



Fig. [6] (k 



(0) 



0.04 and k 



(i) 



0.11): it carries infor- 



mation about the strength of the edge-induced force per 
length, i.e., the chemical edge termination, and the elas- 
tic response of the AGNR which again incorporates the 
effective boundary conditions imposed on the ribbon's 
wavefunctions. 

A striking aspect of Figs. [5] and [5] is the pronounced 
oscillation in N with period three. It translates into the 
cos(27riV/3)-term of Eq. ((2|) that wc take as evidence of 
the boundary conditions feeding back into the ribbons 
elasticity. By contrast, a two-fold periodicit y o riginating 
from even/odd effects is largely suppressed. [43( 

It is a remarkable feature of Eq. ^ that the oscillations 
do not experience significant damping (at least within the 
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FIG. 6. (Color online) Relative deformation of the geometry 
of the unit-cell for hydrogen terminated graphene armchair 
nano-ribbons as a function of the inverse of the width N 1 _ 1 ; 
Values for N = 5 . . . 14 have been considered. The insets high- 
lights the corresponding powerlaw dependencies on the ribbon 
width, which has the remarkable feature that the oscillations 
do not experience significant damping. 



This is in marked contrast to traditional approaches in 
meso- and nanoscopic structures that investigate quan- 
tum size effects near the band-gap by directly probing 
the electronic excitation spectrum, e.g., in transport mea- 
surements. At a length of about 1850 carbon slices an 
N = 12 ribbon will be a full chain of carbon atoms longer 
than an N = 11 ribbon with the same number of chains, 
whereas the N = 13 will be the same amount shorter. 
This difference is large enough that it could be detected 
experimentally. We expect that a termination other than 
hydrogen will only change amount of edge induced stress 
but not the mechanism underlying the elastic response. 
Hence, the values of the constants k will change but not 
the oscillating behavior. 

Financial support by the Center for Functional Nano- 
structures (CFN), and cpu time allocation at the 
OPUS IB and hc3 clusters at the Karlsruhe Institute 
of Technology (KIT) Steinbuch Center for Computing 
(SCC) are gratefully acknowledged. 



system sizes available to our numerics) . In the absence of 
perturbations not included in our model such as ripples 
and other inhomogeneities, we expect that this statement 
should remain valid as long as the deviation between the 
true position of the AGNR-band and the prediction based 
on the selection rules, Eq. ([1]), remains small as compared 
to Ak y . Judging from Fig. UJtliis could be true at least in 
to the regime where the aspect ratio N/L is still small. 

Elastic material responses are usually (approximately) 
volume conserving. It is therefore reassuring to see that 
this is also the case in the present situation. Stretching 
the ribbon in longitudinal direction evokes a transverse 
contraction, see Fig. [SI which eliminates in the leading 
order the strain effect on the volume of the unit cell: 



AA = 



(N-iy 



(3) 



where ka(N) is of the form ^ with constants that can 
be read off the inset of Fig. O n 



(0) 



0.04, k { £ 



0.2. 



Conclusions. We have studied quantum size effects 
in armchair graphene nano-ribbons (AGNR) with hydro- 
gen termination using density functional theory (DFT) 
in Kohn-Sham formulation. By formulating "selection 
rules" that allow to extract (approximately) the elec- 
tronic structure of the AGNR bands starting from the 
four graphene dispersion sheets, we have predicted a 
threefold periodicity of the ribbon's electronic structure 
in the ribbon width N that was confirmed by ab initio 
results. We have also observed how this threefold elec- 
tronic periodicity carries over into the atomic structure 
on AGNR. 

Our results imply that quantum size effects can be 
studied experimentally in H-AGNR by atomic structure 
determination, namely by comparing the length per car- 
bon slice of AGNRs of neighboring widths N, N + 1, . . .. 
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The suppression of period two oscillations is understood 
in the following way. The most important contribution to 
the binding of sp2-hybridized carbon atoms is made by 
in plane a— bonds. Stretching at one edge of the AGNR 
faces a restoring force that originates from the inner 
(j — bonds. Most of this restoring force comes from inner- 
most carbon atoms; only a small contribution is made by 
atoms of the opposing edge. For this reason the elastic 
deformation of the carbon-lattice responding to the H2 
induced stretch on one edge is hardly sensitive to the ter- 
mination of the other edge. 

As long as bond stretching is not too strong, we reside in 
the regime of linear elasticity. Then there is by definition 
no cross-talk between the stimulating stretching forces 
from either edge; the ribbons response is essentially a 
superposition of two independent stimuli, one from each 
edge, each one being largely insensitive to details of ter- 
mination of the other edge. Hence, a two-fold periodicity 
is suppressed. 



